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We investigate the evolution of parton distributions at small x values, ob- 
tained in the case of soft initial conditions. The contributions of twist-two and 
(renormalon-type) higher-twist operators of the Wilson operator product expan- 
sion are taken into account. The results are in very good agreement with deep 
inelastic scattering experimental data from HERA. 

The measurements of the deep-inelastic scattering structure function (SF) 
F2 in HERA El have permitted the access to a very interesting kinematical 
range for testing the theoretical ideas on the behavior of quarks and gluons 
carrying a very low fraction of momentum of the proton, the so-called small x 
region. The reasonable agreement between HERA data and the next-to-leading 
order (NLO) approximation of perturbative QCD that has been observed for 

> 2GeV'^ (seeH and references therein) indicates that perturbative QCD 
could describe the SF evolutiijn up to very low values. 

The standard program to study the behavior of quarks and gluons is 
carried out by comparison of data with the numerical solution of the DGLAP 
equations by fitting the parameters of the x profile of partons at some initial 
Qq. However, if one is interested in analyzing exclusively the small x region 
{x < 0.01), there is the alternative of doing a simpler analysis by using some 
of the exisiing analytical solutions of DGLAP in the small x limit (see, for 
example, uQ). The main ingredients of the studyu are: 

• Both, the gluon and quark singlet densities are presented in terms of 
two components ('-I-' and '— ') which are obtained from the analytical 
dependent expressions of the corresponding ('-I-' and '— ') parton distri- 
butions moments. 

• The '— ' component is constant at small x, whereas the '-f' component 
grows at > Qq as 
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where the LO term a+ = U/jSo and the NLO one b+ = 412//(27/3o)- 
Here the couphng constant as — a^/^ATr), (3q and (3i are the first two 
coefficients of QCD /3-function and / is the number of active flavors. 

We shortly compile below the main results ofi at the leading order (LO) 
approximation and demonstrate some new (preliminary) results, where the 
contributions of higher-twist (HT) operators (i.e. twist-four ones and twist- 
six ones) of the Wilson operator product expansion are taken into account 
in the framework of renormalon model (seel3). The importance of the HT 
contributions at small- a; has been demonstrated inQ. 

1. Our purpose is to show the small x asymptotic form of parton distri- 
butions in the framework of the DGLAP equation starting at some Qq with 
the fiat function: 

f:\Ql) = Aa (hereafter a = q,g), (1) 

where fl^ are the leading-twist (LT) parts of parton (quark and gluon) dis- 
tributions (PD) multiplied by x and Aa are unknown parameters that have 
to be determined from data. Through this work at small x we neglect the 
non-singlet quark component ^. 

The full small x asymptotic results for PD and SF F2 at LO is: 

F2{x,Q^) = e./,(x,Q2), /,(x,Q2) ^ f+ {x,Q^) + f- {x,Q^) , (2) 
where the '+' and '— ' components f^{x, Q^) are given by the sum 

I^M^) = f:''^{x,Q') + ft^^{x,Q') (3) 

of the LT parts /^'^'^{x, Q^) and the HT parts f^^'^ix, Q'^), respectively. 
The small x asymptotic results for PD, f^^'^ 

fp'+{x, Q2) ^ (^A, + io{a) e-~^+(^^> + 0{p) , (4) 

fl^-+{x,Q^) = {[a, + ^aI^pHct) e-^+(i)« + 0{p) , (5) 

where d+(l) = 1 + 20//(27/3o) and d_(l) = 16//(27/3o) are the regular parts 
of d+ and d_ anomalous dimensions, respectively, in the limit n — > 1 ^. The 

" We would like to note that new HERA data! show a rise of F2 structure function at low 
values (Q^ ~ IGeV^) when a; — > (see Fig.l, for example). The rise can be explained in 
a natural way by incorporation of HT terms in our analysis (see Eqs.(^-(^)). 
''From now on, for a quantity k{n) we use the notation k(n) for the singular part when n — > 1 
and k{n) for the corresponding regular part. 
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Figure 1: The structure function F2 as a function of a; for different bins. The experimental 
points are from HI. The inner error bars are statistic while the outer bars represent statistic 
and systimatic errors added in quadrature. The curves are obtained from fits at LO when 
the HT contributions have been incorporated. 

function [v = 0, 1) coincides with the modified Bessel function /i, at s > 
and the Bessel function at s < 0. Using the calculations Ba, we show the 
HT effect in the renormalon case. We present the results only for the terms 
proportional Ojf. some power of ln(l/a;) (full expressions can be found in the 
last papaer ofB), making the following subtitutions in the corresponding LT 
results presented in Eqs.(0)-(0): 
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where ^ and A| ^ are magnitudes of twist-four and twist-six corrections. 
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/;2^-(x,g2) (see Eq.(|)) 
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2. With the help of the above equations we have analyzed F2 HERA data 
at small x from the HI collaboration. We have fixed the number of active fla- 
vors /=4 and Aj^(n/ = 4) = 250 MeV, which is a reasonable value extracted 
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from the traditipiial (higher x) experiments. Moreover, we put Ai a = A2,a in 
agreement withli3. 

The results are shown on Fig. 1. We found very good agreement between 
our approach based on QCD and HERA data. The (renormalon-type) HT 
terms lead to the natural explanation of the rise of F2 structure function at 
low values of and x. 

As a next step of our investigations, we plan to finish this study and to 
investigate HT contributions to PD and SF relations, observed inEJIlj. 
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